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Abstract 

o 

^NJ ■ We study the reaction-diffusion system, its stationary 

O . solutions, the behavior of the system near them and dis- 

D | cuss similarities and differences for different boundary con- 

ditions. 

m 
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< 

ji ■ 1 Introduction 

We deal with a nonlinear parabolic system of the form 
(1.1) u t = DAu + /(«), 

> 

\Q ; with Neumann homogeneous boundary condition 

d.2) £i„ = o 

or, incidentally, Dirichlet homogeneous boundary condition 
(1-3) «|an = 0, 

in a bounded domain f2 C IR m with boundary <9f2 being m — 1- 

dimensional sufficiently smooth manifold. Here, we consider only 
c3 | 
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classical solutions u : O — > R™ (we emphasize that we have vector- 
valued functions, since ( II .ip is, in fact, a system of n parabolic 
equations. We denote by D := diag(c? 1 , . . . , d n ) a diagonal matrix 
with positive diagonal entries di, i = 1, . . . , n, and by / : M™ — > R n 
a continuous function defined on the cone of nonnegative vectors 
x = (xi, . . . , x n ) G W 1 , where Xi > 0, i — 1, . . . , n. Such parabolic 
systems are largely investigated since they model kinetics of chem- 
ical reactions - each coordinate of u measures density one of inter- 
acting components, or they model development of several biological 
species living on the same area Q and interacting in different ways 
(preys and predators, symbiosis, competition). We keep in mind the 
second application that causes Neumann condition is more natural 
(comp. [2]). System (II .ltil.2|) defines a semiflow $ t , t > 0, on an 
appropriate space X a , if one uses the theory of sectorial operators 
(comp. [61 [9]) or a semiflow on a cone of nonnegative continuous 
functions, if one applies the theory of monotone dynamical systems 
(comp. [IE]). We work in spaces of continuous functions, since our 
main example of / = (/i, . . . , /„) is 

(1.4) fi(u) = m ^1 - y~] ajjU-^j , 

where all coefficients are positive. If / is a C 1 -function such 
that 

fi(u) = for Ui = 0, i = l,...,n, 

and 

(1.5) f£<0, for i^j, 

then (11.11) defines a semiflow $ t ? t > 0, on the cone of nonneg- 
ative continuous functions u : Q — > [0, oo) n which is competi- 
tive in the sense of Hirsch (see [T5| [T^]). In particular, it means 
that all solutions of the system starting with nonnegative functions 
ip = u(-,0) > are global in time and nonnegative for any t > 0. 
Moreover, if 

U(;T)<U(;T) 

are two such solutions comparable at any time T > 0, then 
u(-,t) < u(-,t) for t < T. 
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2 Steady-state solutions, single species 



First, we are interested in steady-state solutions, i.e. time indepen- 
dent solutions. They satisfy elliptic system: 

(Jll 

(2.1) Au = -D-'fiu), — \ dn = . 

The obvious examples of such functions are zeros of the nonlin- 
ear term /. An important question is if there are not exist other 
steady-state solutions. They correspond nonuniform distributions 
of populations on the environment Q which do not change in time. 
We can exactly investigate the case of one species (n — 1). If Q C M., 
the analysis is standard and rather simple. 
(2.2) 

u" = -D^uil - u) in (0,L), u'(0) = = u'(L), u>0. 

The second order ODE is an example of a conservative system with 
one degree of freedom so it has a first integral (energy) 

E(u,u') = 1 . 

v ' ; 2 2D 3D 

It enables to find the phase portrait of the system: 



u' 



One can see that there is no nonnegative solutions for any dif- 
fusion coefficient D and any length L of the environment. On the 
other hand, if one use the Dirichlet boundary condition u(0) = = 
u(L), then a priori such solutions can exist. They will correspond 
trajectories cutting axe u' twice. One can compute the time (we 
use standard notions from the theory of dynamical systems how- 
ever, here, the independent variable is interpreted as a spatial one) 
between two passes of this axe. Denote (after [T2]) by F the real 
function F(u) = u 2 /2 — u 3 /3 and let \i be the value of u where the 
trajectory cut axe u. Then this time equals 

du f 1 jidz 



Jo JF u) - F u) Jo 



y/F(fJ.) - F(U) Jo y/F{fJ,)-F(fJ,z)' 

It is obvious that \x changes between and 1, that L is an increasing 
function of \x tending to oo as \i — > 1~. One can also compute the 
limit 

f 1 dz 

lim L(u) = 2VD / ; = irvD. 

It can be interpreted that there is exactly one nonconstant steady- 
state solution iff the length of the environment L is greater than this 
limit Lkiss — called the KISS size. If L < tt\/D, then there 

are only constant steady-state solutions u = and u = 1. From 
another point of view, if the length L is fixed, then nonconstant 
time-independent solutions exist when the diffusion coefficient D is 
sufficiently small. 

Now, we study the case, when Q C M m with m > 1. The most 
interesting dimension from the biological point of view is m = 2. 
Q = B(0,R) (the disk centered at with radius R) is the simplest 
set and we can easily look for radial solutions of (12. ip . Our analysis 
can be repeated in larger dimension without troubles. If we denote 
by v! the derivative with respect to the radial coordinate, then we 
get the following boundary value problem: 

(2.3) u" + -u' = - u) in (0, R), u'(0) = = u\R). 

For Dirichlet's boundary condition u\dQ = 0, we have for radial so- 
lutions: u'(0) = = u(R). For both problems, we look at solutions 
of the second order ODE with initial values u(0) = c, u'(0) = 0, 
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where c > 0. It is easy to see that for c > 1 the right hand side of 
ODE (12. 3p is positive, thus the solution increase to infinity in finite 
time and cannot satisfy the Dirichlet nor Neumann conditions at 
R. If c G (0, 1), then the solution is concave near 0, u decreases but, 
by comparison with equation f l2.2|) . slowler then in the last system. 
It means that, for any c, u{r\) — when r\ > Lkiss/2- Thus, for 
Dirichlet 's boundary condition, we have the same phenomenon as 
in dimension 1, but with larger KISS size. Similarly, there is no 
nonconstant solutions for Neumann's problem. Below, we present 
numerical solutions for initial problems 

u" H — u = ——u(l — it), u(0) = c, u'(0) = 
r D 

with several c > and D = 0.1. 




3 Steady-state solutions, two species 

For two species the situation is completely different. Consider the 

system 

(3.1) 

f u t = d 1 Au + u(l-u-bv), pU = 0, u{-,0) = (p 

5 



where interspecies-competition coefficients b, c are positive constants, 
diffusion coefficients d±,d 2 are positive and Q is an open bounded 
set in IR m with smooth boundary. We study only nonnegative so- 
lutions. There are always three equilibria constant in space and 
time: 



Volterra ODE with the above equilibria. The dynamical system 
given by this ODE is easily investigated: Pq is its repeller, P u is an 
attractor if c > 1, P v is an attractor if h > 1, Pi is global attractor 
(for nonnegative solutions) if b, c < 1. More exactly, if both b and c 
are greater than 1, then Pi is a saddle point and its stable manifold 
W S (P\) is the sum of a heteroclinic trajectory from P to Pi and a 
trajectory from infinity to Pi. This manifold cuts the set u, v > in 
two sets: all trajectories starting from the set containing P u tend to 
P u , trajectories starting from the second set tend to P v . The method 
of monotone dynamical systems enables us to state a similar result 
for full parabolic system (13. ip : 

(i) if b, c < 1, then all solutions with nontrivial <p, ip > tend to 
Pi as t — t- oo; 

(ii) if b < 1 and c > 1, then all solutions with nontrivial ip, ip > 
tend to P u as t — > oo; 

(iii) if b > 1 and c < 1, then all solutions with nontrivial <p, ip > 
tend to P v as t — > oo; 

(iv) if b, c > 1, then the stable manifold W S {P{) is again one- 
dimensional, P u attracts solutions with ip > Uo, ip < fo for some 
(uq, v o) G W s (Pi), and P v attracts solutions with p < uo, ip > v 
for some (uo,Vo) G W S (P\). 

The proof can be found in [TH] . 

The last case is the most interesting since we do not know the 
dynamics for any (</?, ip) not comparable in the above sense with any 
function from W^Pi). For example, some nonconstant equilibria 
can exist and even they can be asymptotically stable. However, we 
can use the following result due to Conway, Hoff and Smoller [3] in 
many cases: 



P = (0,0), P U = (1,0), P„ = (( 
and if b, c < 1 or b, c > 1, the fourth equilibrium 



(0,1) 




If d\ — di — 0, there is no diffusion and we have the standard Lotka- 
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Let u be a solution to 



u t = DAu + f(u), —\ dn = 0, u(.,0) = uo, 

where u = (u±, . . . , Mm); -D is a symmetric positively definite matrix, 
Q is an open bounded subset of R n with smooth boundary, / is a 
C 2 function and uq G L 2 (fi). Assume that there exists a bounded 
positively invariant set E C K m such that -u takes values in E. 
Denote by A the first eigenvalue of —A with |^ \ dn = (here, one 
can take m = 1, and we mean the first positive eigenvalue - is 
also eigenvalue), by d the lowest eigenvalue of matrix D, M : = 
sup ueE and at last 

a := Xd - M. 

Then there are four positive constants q, % — 1, 2, 3, 4 such that 
(3.2) ||V«(-,t)IU»(n) < Cl e- CTt , 



(3.3) I !«(•»*) - "(Oilmen) < c 2 e' 



where m is a spatial average of u and it satisfies vl = f(u) + g(t), 
u(0) = f n u / /Lt(Q) , g a is a function satisfying |g(£)| < c 3 e" CT *, 



(3.4) | \u(-, t) - u(i) | | L oo ( n) < c 4 e 



if D is a diagonal matrix. 

In our investigations, E = {(u, v) : -u, i> > 0, tt < 1 — bv or v < 
1 — cu} if (ip, ip) takes values in this set or E is the smallest tri- 
angle with vortices P , (a, 0), (0, a) containing (sup if, sup ip). Since 
the first set attracts all solutions, it contains all equilibria. If one 
has a > 0, then lim^oo u(t) = const and we the solution u(x, t) 
tends to this constant as t — > oo. It follows that there is no noncon- 
stant in space steady-state solution. After easy though laborious 
computations one can find the constant M. We have 



f'(u,v) = 



1 — 2u — bv —bu 

—cv 1 — 2v — cu 



\\f'(u,v)\\ 2 = {b 2 +c 2 +4){u 2 W)+4:{b+c)uv-4{u+v)-2(cu+bv)+2 

and the maximum of the last function on the set E equals 2 (it is 
reached at the origin). Thus, we have obtained 
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Theorem 1. If both diffusion constants di,d,2 are sufficiently large, 
namely 

min((<i, d2) > — — 
A 

where A is the first positive eigenvalue for —A with Neumann ho- 
mogeneous condition, then there is no nonconstant steady-state so- 
lution. For the case Q = (0, L) C M. as in the previous section, we 

2 

have X = j2 and we have no nonconstant equilibrium if 

L < 2~ 1/4 7rmin(di,<i2). 

Compare this number with the KISS size from the previous sec- 
tion and notice that, here, a priori we have nonconstant equilibria 
for the Neumann boundary condition. There is a numerically stud- 
ied example of Matano and Mimura [16] where Q is a set in the 
plane consisting two squares joined by a thin strip and b, c > 1 
such that there is a nonconstant positive equilibrium. Moreover, it 
is asymptotically stable. If domain Q is convex all equilibria are 
stable (see [H]), hence it is not surprising that in this example the 
domain is such. On the other hand, one can understand ecological 
sense of the shape of the domain: in the first square the first species 
wins, in the other the second one; the thin strip makes migrations 
between squares more difficult hence both species can coexist. 

Most results from the theory of monotone dynamical systems 
one used concern systems in ordered Banach spaces called SOP in 
the monograph by H. Smith [18]. Competitive parabolic systems 
are SOP if we use order: 

(u,v) -< (u,v) u(x) < u(x), v(x) > v(x) for any x. 

The semiflow generated by (13. ip preserves this order, i.e. if (ip, ip) -< 
(<f, tp), then t),v(-,t)) -< (u(-,t),v(-, t)) for any t > 0. It enables 
us to get some information on w-limit sets of our system. This 
choice of the order is typical for competition of two species and can 
be explain qualitatively by ecological arguments: if the first species 
dominates the second one and one increases the population of the 
first species and decreases of the second one then the domination 
conserves. 
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4 Steady-state solutions, many species 



For more than two species the situation is much more complicated. 
If n = 3 then, roughly speaking the first species can dominate the 
second one, the second one can dominate the third one and this 
last species can dominate the first one. The simplest mathematical 
model of this case is given by May and Leonard [T7] for ODE: 

{x = x(l — x — ay — (3z) 
y = y(l-Px-y- az) 
z = z(l — ax — f3y — z) 

with a, (3 > 0, a + (3 = 2. The typical w-limit set is a limit cycle 
but there are three heteroclinic trajectories joining three equilibria 
(1, 0, 0), (0, 1, 0) and (0, 0, 1). We have studied a slightly more gen- 
eral case in [TU] however the behavior of the system is very similar. 
Consider the parabolic system for three species: 

U t = DAU + f(U) } 
(4-2) f \an = 0, 

U(;Q) = tp, 

where U = (u,v,w), D = diag((ii, o?2, ^3), 

u(l — diu — b\v — C\w) 
v(l — a 2 u — b 2 v — c 2 u>) , 
w(l — a 3 u — b 3 v — c 3 w) 

all constants in the above formulas di, cii, bi, Ci are positive, x G 
tt C R m . 

We study the system (14.21) and its spatially homogeneous ODE 
system 

(4.4) U' = f(U) 

in the open set 

V := {(u, v, w) : u,v,w > 0} 

which is obviously invariant and the same is true for its closure T>. 
Fixed points of the system (I4.4p can be easily found - four of them 
always lie in D : 

P = (0, 0, 0), P u = (or 1 , 0, 0), P v = (0, b 2 \ 0), P w = (0, 0, C3 1 ), 



(4-3) f(U) = 
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next four are 

b 2 - h 



P 



cii — a 2 







aib 2 - a 2 bi ' aib 2 - a 2 bi 

c 3 - c 2 



P„ 



C 3 — C\ 0\ 



p 

J. ny 



0. 



(IIC3 — a 3 c\ 
b 2 - b 3 



aic 3 



b 2 c 3 - b 3 c 2 b 2 c 3 - b 3 c 2 
and Pi = (a, (3, 7) being the solution of the equation 





u 




' 1 " 




cii 




C\ 


M 


V 




1 


where M = 


a 2 


b 2 


c 2 




w 




1 




a 3 


h 


c 3 



We assume that W := det M^fl what means that none of two of 
planes 

Hi : diU + biV + CiW = 1, i — 1, 2, 3 

are parallel. It is easy to see that f'(Po) = I and its unique eigen- 
value 1 is positive, thus P is a source for (14.41) . Similarly 



f'(Pu) 



— 1 —b\ja\ —ci/ai 
1 — a 2 /ai 
1 — 03/01 



and if max(a2,as) < a l5 then one eigenvalue —1 is negative and 
two remaining ones are positive. The corresponding eigenspaces 
are: {(w, 0,0) : u G M} and {(0,v,w) : v,w G R}. The stable 
manifold of P u is 1-dimensional: 

{(m,0,0) : u > 0} 

and by the Stable and Unstable Manifold Theorem all trajectories 
of (14.41) except starting in the stable manifold cannot tend to P u as 
t — > +00. Similar arguments work for P v and P w under conditions: 

max(6i,6 3 ) < b 2l max(ci,c 2 ) < c 3 . 

In [TU], we proved even more: 

Lemma 1. // 

(4.5) min(a 2 ,a 3 ) < «i, min(6i,6 3 ) < b 2} min(ci,c 2 ) < c 3 , 

then P u , P v and P w do not belong to the u-limit set uj{P) of any 
point P G P. 
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Dividing the set T> into three pieces: 

V + := {(u, v, w) G T> : min(ajM + b{V + Ciiu) > 1}, 

i 

A := {(u,v,w) G V : min(ajM+6jf +Ciw) < 1 < max(ajti+6jf +c i w)} ) 
T?_ := {(u, v , w) G V : max(ajU + Ojt> + qw) < 1} 

i 

(X> + (resp. P_) is the set of points sitting under (resp. over) all 
three planes H i} i = 1, 2, 3, A = V\ (X? + U£>_)) we got \W\ following 
result for ( Oj) : 

Lemma 2. T/ie set A zs positively invariant and all trajectories 
in T> eventually come into A. Moreover, A contains any compact 
invariant set that contains no fixed points. 

Thus, u(P) C A for any P e V. Notice that P u , P v , P w G A 
and similarly P uv , P uw , P vw if they belong to V. 

The last fixed point Pi = (a, (3, 7) can lie in T> (and then in A) 
or outside this cone. The following theorem excludes the existence 
of a periodic trajectory for ( 14. 4p in T> if Pi G T>. 

Lemma 3. (fTE/, p. Prop. 4-3) Let T be a nontrivial periodic 
trajectory of a competitive system in T> C M 3 and 

rcb,?] := {£: Pi<&<fc, i = l,2,3} CD. 

TTien t/ie set 0/ a// points x which are not related to any point 
y G T (relation x < y means < Wj /or any i j has two components, 
one of them is bounded and contains a fixed point. 

Hence, if we want to have a nontrivial periodic trajectory, then 
Pi must belong to T> and we have two options: 
(i) W > and three other determinants 







" 1 


61 


Cl 




:= det 


1 


b 2 


c 2 






1 


h 








di 


1 


Cl 


w v 


= det 


a 2 


1 


c 2 






% 


1 


C 3 
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W w := det 



ai bi 1 
a 2 b 2 1 
a 3 6 3 1 




are positive or 

(ii) W < and the above three determinants are negative. 
The main result of [10J is the following 

Theorem 2. Assume Let all four determinants W, W u , W v , 

W w be positive and p := a^a + b 2 (3 + C37 — 1 < 0. Then, for any 
point P G T> that does not belong to the half-line 



s > 0, 



the ui -limit set u(P) for ft4-4\ ) ^ s a periodic trajectory. For P from 
this half-line, u (P) = P\. 

If we combine this theorem with the result [3] cited in the pre- 
vious section and a classical result of L. Markus fo] that cu-limit 
sets of (14.21) are the same as corresponding (14.41) . we obtain 

Theorem 3. Suppose that inequalities ( |^.<5| ) hold, four determi- 
nants W, W u , W v , W w are positive, p < and all diffusion coeffi- 
cients di are sufficiently large, then for any continuous, nonnegative 
function (p with all coordinates nontrivial, either 

lim U(x, t) = Px 

t— >oo 

or there exists a periodic function U : R — > A C V such that 

lim |17(2;, *) - U(t)\ = 

t— >oo 

and in both cases limits are uniform with respect to x. The first 
case takes place only if j n <p belongs to the half-line starting from 
the origin. 

The condition for the diffusion coefficients is given by: 
Aiming: i = 1, 2, 3} > sup{|| /'([/) || : U G A} 
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where Ai is the first positive eigenvalue of —A with Neumann's 
boundary condition, but numerical experiments show that the as- 
sertion of the above theorem holds also for essentially larger these 
coefficients and even if only one of them is great and two others 
are very small. However, for all coefficient being small there ex- 
ists spatially heterogeneous steady-state solutions which attracts 
other solutions. Obviously, these functions take values in positively 
invariant subset A. 

If the number p defined in Theorem [2] is positive, then all eigen- 
values of f'(Pi) have negative real parts - calculations in [10J show 
that the characteristic polynomial equals 

Q(X) = (A + 1)(A 2 + p\ + apryW). 

It follows that 

lim U(x, t) = P 1 

t— >oo 

uniformly in x G Q for any starting point ip as in the last theorem. 
The case p = is difficult for investigations since the behavior of 
the system near Pi cannot be found by the linearization; Pi is not 
hyperbolic. In the special case from the paper of May and Leonard 
[TT] ai = &2 = C3 = 1, 61 = c 2 = a 3 , ci = a 2 = 63 if p = there are 
two first integrals of (14.41) and trajectories can be found apparently. 
The whole triangle spanned by P u , P v and P w is fulfilled by limit 
cycles around Pi. For sufficiently small diffusion coefficients these 
limit cycles describe the asymptotic behavior of the system (14. 2p . 

If Pi does not belong to V, then all solutions tends to one of 
the others steady-state solutions at least for small diffusion coeffi- 
cients, when there is no nonconstant steady-state solutions due to 
the above mentioned arguments. It means that at least one of the 
species extincts. 

5 Some numerical simulations 

The most interesting case considered in the previous section is pre- 
sented in the theorem: almost all solutions of (j4.2p tend to periodic 
functions given by the system of ODEs (14.41) . We can investigate 
numerically equation (14. 2 j) with Q = (0, 1) C R and / of the form 
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from the previous section. Put matrix M of competition coefficients 



M 



2 1.1 
3.1 2 
0.95 2.9 



3.1 
0.9 
2 



This choice ensure the assumptions of Theorem 3 hold. The first 
positive eigenvalue of this degenerate Laplacian equals Ai = tt 2 and 
one can compute the maximal value of the norm of the derivative 
\\f'(U)\ \ on the set A : 

\\f(U)\\<V3. 



Thus the critical value of diffusion coefficients is d := For 
mm{di : i — 1,2,3} > d, almost all trajectories tend to periodic, 
spatially constant functions. Since the theorem of Conway, Hoff 
and Smoller gives only the sufficient condition, it is not surprising 
that even for smaller di's the assertion is true. 

We have found numerically (by using Maple 10) solutions of 
(14.21) with the above matrix M and 



di = 10 



-3 



do = 2 • 10 



-3 



d* = 0.5 • 10" 



<p(x) = [6x 2 (l 



x 



x ; 



\2x\l 



x 



The type of this initial function is natural if we are seeking classical 
solutions: the normal derivative of the initial function should vanish 
at both boundary points x — 0, x = 1. Below, we present the plot 
of the graph of the second coordinate of U for three values of x : 
0.1, 0.5 and 0.9 as the function of time t for the range [0, 100] : 
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u 20 40 60 80 100 

t 

The plots for three values of x are completely different, hence 
the solutions of the parabolic system tend as t — > oo to functions 
which are not spatially constant. Nevertheless, they seem to be 
periodic as functions of time. Compare this plot with the plot of 
the second coordinate v of the solution to ODE (14.4ft with initial 
point (0.1,0.0095238,0.0333333) which is the spatial average of <p. 
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20 40 60 80 100 



6 Stability of steady-state solutions 

The stability analysis in the linear approximation for constant solu- 
tions of the parabolic system is standard (see, for example [20], 19]). 
Consider (14. 2 1) not necessarily with U taking values in M 3 but in 
R n . If P is a zero of /, then U = P is a solution both 04. 2 p and 
(I4.4I) and it is asymptotically stable for both systems again if all 
eigenvalues of f'(P) have negative real parts. Below, we shall study 
the analogous problem for periodic solutions of f)4.4p . 

Let t \-t U(t) be such a solution. Denote by A(t) = f(U(t)), 
by Afc, k = 0, 1, 2, ... , the sequence of all eigenvalues of —A with 
Neumann boundary conditions (A = < Ai < A 2 < . . .) and 
by x i — y efc(x) corresponding eigenf unctions. Since the spectrum is 
discrete and the operator is self-adjoint, the eigenf unctions form a 
complete orthonormal system in L 2 (Q). One can solve the linearized 
problem 

(6.1) U t = DAU + A(t)U, ^U = 0, U(;Q) = <p 
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by using the Fourier method. If we put 



oo 



U(x,t) = y]e fc (x)g fc (t) 



fc=0 



in (16. ip . we have Neumann's condition satisfied and vector- valued 
function g% should be a solution of the initial problem: 



for any k, where is a coefficient of the Fourier expansion of ip 
with respect to the orthonormal system {e^ : k = 0,1, . . .}, i.e. 



Since matrix D is diagonal, it commutes with A(t) and multipliers 
(Floquet theory) of — XkD + A(-) are of the form exp(— Xkd)g, where 
d belongs to the interval [min{dj}, max{<ij}] and g is a multiplier 
for matrix A(-). Hence decays exponentially as t tends to +00, if 
\g\ < exp(Afc min{<ij}) for any multiplier g. It is well known that one 
of multipliers for A equals 1, thus this inequality cannot hold for 
k = as Ao = 0. But, in spite of this, if the remaining multipliers of 
A sit in the open unit disc and 1 is simple, then solution U of (14.21) 
is orbitally asymptotically stable - see [20] Theorem 8.2.3, p. 251, 
i.e. there exists a neighborhood of the periodic orbit T := {U(t) : 
t G M} such that solutions U\ with initial function ip taking values 
in this nhbd tend to U in the sense 



where the above norm means usual one in H 1 and to is a number de- 
pending on ip called asymptotic phase. For 1-dimensional domains 
Q, one has the canonical embedding of H 1 and the above limit is 
uniform in x. For more natural environments (Q C M 2 ), we cannot 
use this argument to get the uniform limit. 

7 Concluding remarks 

Recently, a lot of important papers on spatial heterogeneity models 
of competiting species [H El El El d3l E] , for instance. They consider 



(6.2) 



9 , k(t) = (-hD + A(t))g k {t) 



9k(0) = c k 




lim \\U 1 (;t)-U(t + t )\\ = 0, 



17 



mutual interplay of diffusion and competition which gives many 
interesting phenomena for such systems. However, most of results 
are obtained for two species but spatial heterogeneity is included 
in the model - some coefficients depend on variable x. It seems 
that some new effects can be obtained if we study interaction of 
more than two species. The present author hope that some possible 
directions of such investigations have been indicated above. 
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